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Abstract 

We investigate strong-to-weak coupling transitions in D = 2 + 1 SU{N —>■ oo) gauge 
theories, by simulating the corresponding lattice theories with a Wilson plaquette action. We 
find that there is a strong-to-weak coupling cross-over in the lattice theory that appears to 
become a third-order phase transition at = oo, in a manner that is essentially identical to 
the Gross- Witten transition in the D = 1 + 1 SU{oo) lattice gauge theory. There is evidence 
of an additional second order transition developing at iV = oo at approximately the same 
coupling, which is connected with Z^v monopoles (instantons), thus making it an analogue of 
the first order bulk transition that occurs in D = 3 + 1 lattice gauge theories for > 5. We 
show that as the lattice spacing is reduced, the N = oo gauge theory on a finite 3-torus suffers 
a sequence of (apparently) first-order Zn symmetry breaking transitions associated with each 
of the tori (ordered by size). We discuss how these transitions can be understood in terms of 
a sequence of deconfining transitions on ever-more dimensionally reduced gauge theories. We 
investigate whether the trace of the Wilson loop has a non-analyticity in the coupling at some 
critical area, but find no evidence for this. However we do find that, just as one can prove 
occurs in D = 1 + 1, the eigenvalue density of a Wilson loop forms a gap at = oo at a 
critical value of its trace. The physical implications of this subtle non-analyticity are unclear. 
This gap formation is in fact a special case of a remarkable similarity between the eigenvalue 
spectra of Wilson loops in = 1 + 1 and D = 2 + 1 (and indeed D = 3 + 1): for the same 
value of the trace, the eigenvalue spectra are nearly identical. This holds for finite as well as 
infinite A^; irrespective of the Wilson loop size in lattice units; and for Polyakov as well as 
Wilson loops. 



1 Introduction 



A phase transition is associated with a singularity in the partition function, and so requires 
an infinite number of degrees of freedom. Usually that requires an infinite volume. One of the 
peculiarities of large- field theories is that one can have phase transitions on finite, or even 
infinitesimal, volumes at = cxd because in this case we have an infinite number of degrees 
of freedom at each point in space. The classic example in the context of gauge field theories 
is the Gross- Witten transition P that occurs in the D = 1 + 1 SU (oo) lattice gauge theory 
(with the standard Wilson action). In this case the theory is analytically soluble and one finds 
a third order phase transition at A^ = oo |T] at a value of the bare coupling that separates 
the strong and weak coupling regions. The theoretical and practical interest of such phase 
transitions, particularly in D = 3 + 1, has recently been reviewed in 0. 

In D = 3 + 1 SU (N) gauge theories numerical studies reveal the existence for A^ > 5 of a 
first order 'bulk' transition separating the weak and strong coupling regions [SHHEj- One also 
finds that the deconfinement transition, which is first order for A^ > 3, becomes sharper on 
smaller volumes as A^ increases suggesting [0] that here too one will have a phase transition on 
a finite volume at A^ = cxd. Indeed there appears to be a whole hierarchy of finite volume phase 
transitions at A^ = cxd [71 12] which are, we shall argue below, related to the deconfinement 
transition. 

These are all in some sense strong to weak coupling transitions, and this has led to the 
conjecture [SI 121 that Wilson loops in general will show such N = oo transitions as the 
lattice spacing decreases, when the physical size of the loop passes some critical value. Such 
a transition in D = 3 + 1 could have interesting implications for dual string approaches to 
large-A^ gauge theories, as well as providing a natural explanation for the rapid crossover 
between perturbative and non-perturbative physics that is observed in the strong interactions 
m El In fact it is known [HI CHI that in the N = ooD = l + l continuum theory the 
eigenvalue spectrum of a Wilson loop suffers a non-analyticity for a critical area that is very 
similar to that of the plaquette at the Gross- Witten transition. However, in contrast to the 
Gross- Witten transition, there is no accompanying non-analyticity in the trace of the loop 
and it is unclear what, if any, are its physical implications. 

In this paper we investigate the existence of such phase transitions in D = 2 -|- 1 SU{N) 
gauge theories, as a step towards a unified understanding of these phenomena in all dimensions. 

In the next Section we briefiy describe the SU (N) lattice gauge theory and how we simu- 
late it. There follows a longer section in which we review in more detail what is known about 
the large- A^ transitions and, in some cases, we extend the analysis. (We are interested in tran- 
sitions that may be cross-overs or actual phase transitions, and when we refer to 'transitions' 
in this paper it may be either one of these.) Having established the background, we move on 
to our detailed numerical results. Our conclusions contain a summary of our main results. 
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2 SU(A^) gauge theory on the lattice 



We discretise Euclidean space-time to a periodic cubic lattice with lattice spacing a and size 
LqX Li X L2 in lattice units. We assign SU{N) matrices, Ui, to the links / of the lattice. (We 
sometimes write Ui as U^{n) where the link I emanates in the positive /i direction from the 
site n.) We use the standard Wilson plaquette action 



5 = /35^{l-lReTr[/,} 



N 

V 

where f/p is the ordered product of the S\J (N) link matrices around the boundary of the 
plaquette p. The partition function is 

2N 

r\dUiexp{-S) ; lim/? = — . (2) 

Exactly the same expression defines the lattice gauge theory in D = 1 + 1 and D = 3 + 1 except 
that P = 2N/a?g'^ and P = 2N/g^ respectively. Eqn© also defines the finite temperature 
partition function, if we choose 

T=^ ; Li,L2>Lo. (3) 

We simulate the above lattice theory using a conventional mixture of heat bath and over- 
relaxation steps applied to the SU{2) subgroups of the SU{N) link matrices. 

It will sometimes be convenient to distinguish couplings, inverse bare couplings and (crit- 
ical) temperatures in different space-time dimensions, D, and we do so using subscripts or 
superscripts, e.g. gjj, Pd, . Where there is no ambiguity we will often omit such sub- 
scripts. 

To obtain a smooth large N limit we keep g'^N fixed fT^. It is therefore useful to define 
the bare 't Hooft coupling. A, and the inverse bare 't Hooft coupling, 7, 

A = «/iV , 7 = i = Jj. (4) 

Various numerical calculations have confirmed that a smooth N 00 limit is indeed obtained 
by keeping A fixed, both in D = 2 + 1 ^Sj and inD = 3 + l IHIHl ^3 ^^"^ ^^^^ keep the 
cut-off a fixed as — > cxo one should keep 7 fixed. 

A useful order parameter for finite volume phase transitions is provided by taking the 
Polyakov loop, which is the ordered product of link matrices around the /x-torus, and 
averaging it over the space-time volume: 

= Cm XI iv'^^ I n ^t^(^o, nu ^2) > (5) 

where the normalisation is = W^^^Ly. When the system develops a non-zero value for 
(T^) this indicates the spontaneous breaking of a global Zn symmetry associated with the 
yU-torus. In particular such a symmetry breaking occurs at the deconfining temperature, if the 
yU-torus defines the temperature T. 
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3 Background 



3.1 The ' Gross- Witten' transition 

By fixing gauge and making a change of variables, one can sliow [1 tliat tlie partition function 
of the D = 1 + 1 SU{N) lattice gauge theory (with the Wilson plaquette action) factorises into 
a product of integrals over SU{N) matrices on the links and the theory can be explicitly solved. 
One then finds a cross-over between weak and strong coupling that sharpens with increasing 
N into a third order phase transition at = oo. In terms of the plaquette, Up = KeTrUp/N, 
this shows up in a change of functional behaviour 



More detailed information about the behaviour of plaquettes and Wilson loops can be gained 
by considering not just their traces but their eigenvalues. The eigenvalues of an SU{N) 
matrix are just phases, A = expjio;}, and are gauge-invariant. (We also use A for the 't Hooft 
coupling: which is intended should be clear from the context.) As /? ^ the eigenvalue 
distribution p{a) of a Wilson loop becomes uniform while as /? — oo it becomes increasingly 
peaked around a = 0. As shown in pp, at the Gross-Witten transition a gap opens in the 
density of plaquette eigenvalues: in the strongly-coupled phase the eigenvalue density is non- 
zero for all angles — tt < a < tc, but in the weakly-coupled phase it is only non-zero in the 
range —ac < a < etc, where < vr [J. 

In D = 3 + 1 it is known that at A^ = cxo and indeed for A^ > 5 jH E], there is 
a strong first order transition as f3 is varied from strong to weak coupling. Calculations 
in progress suggest that the plaquette eigenvalue distribution does indeed show a gap 
formation at A^ = oo that is similar to the D = 1 + 1 Gross-Witten transition. However the 
first order transition itself is usually believed to be a manifestation of the phase structure 
one finds with a mixed adjoint-fundamental action as discussed below. This finite-A^ phase 
transition 'conceals' any underlying N = oo Gross-Witten transition and makes the latter 
hard to identify unambiguously. 

In D = 2 + 1 there has been, as far as we are aware, no systematic search for a Gross-Witten 
or 'bulk' transition, and this is one of the gaps that the present work intends to fill. 

3.2 Wilson loop transitions 

The Gross-Witten transition involves the smallest possible Wilson loop, the plaquette. On the 
weak coupling side the plaquette can be calculated in terms of usual weak-coupling pertur- 
bation theory; but this breaks down abruptly at the Gross-Witten transition, beyond which 
a strong coupling expansion becomes appropriate Ij. The coupling is the bare coupling and 
hence a coupling on the length scale of the plaquette. Thus one might interpret the transition 
as saying that as one increases the length scale, there is a critical scale at which perturbation 
theory in the running coupling will suddenly break down. 




A > 2 



A < 2. 



(6) 
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One might imagine that this generahses to other Wilson loops: i.e. when we scale up 
a Wilson loop, at some critical size, in 'physical units', there is a non-analyticity. In fact 
precisely such a scenario has been conjectured for SU{N oo) gauge theories in D = 3 + 1 
|Hl 12] • Unlike the lattice Gross- Witten transition, this would be a property of the continuum 
theory. 

Such a non-analyticity does in fact occur for the SU {N oo) continuum theory in D = 
1 + 1 CD] . The transition occurs at a fixed physical area 

Very much larger Wilson loops have a flat eigenvalue spectrum p{a) which becomes peaked 
as A ^crit- A decreases through Acrit a gap appears in the spectrum near the extreme 
phases a = ivr. So for loops with A < A^.^^ the eigenvalue density is only non-zero for 
— «c where etc < vr, and etc as A 0. The non-analyticity at A = Acrit is in 

fact more singular than for the Gross- Witten transition, in that the derivative dp/ da diverges 
at etc = ivr. However, unlike the Gross- Witten transition this is not a phase transition: the 
partition function is analytic. Moreover the trace of the Wilson loop, and the traces of all 
powers of the Wilson loop, remain analytic in the coupling. Thus it is unclear what if any is 
the physical significance of this non-analyticity. 

In this paper we shall investigate whether such a non-analyticity develops in = 2 + 1 
SU (N) gauge theories and whether it is accompanied by any non-analyticity of the trace. The 
implications could be very interesting ^T] and this makes a search in D = 2 -|- 1 (and even 
more so in D = 3 -|- 1 fE]) well worth while. 

3.3 Mixed actions 

In D = 3 -|- 1 the strong-to-weak coupling transition occurs already at finite A^. It is a cross- 
over for < 4 and is first order for > 5 [H E] • The conventional interpretation of this 
'bulk' transition proceeds by considering a generalised lattice action containing pieces in both 
fundamental and adjoint representations [T^ : 

S = /?/E{l-^I^^Tr,f/,|+/3;J]|l-^^Tr.f/,| 
= E {l - J^^^T^'fUp] + E {l - ]^Tr;f/tTr^^,| 



(8) 



where we have used Tr^f/p = |Trjf/pp — 1. By considering smooth fields one finds that at 
weak coupling, and to leading order in g^, one obtains constant physics by keeping constant 
the linear combination /?/ + 2f3a- 

Consider the limit (3a ^ oo while keeping /?/ fixed. This requires \Ty fUpl"^ / N"^ = 1 which 
implies that the link matrices are elements of the centre. Fluctuations between different 
elements of the centre are controlled by the linear plaquette term multiplied by /?/ which 
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means that what we have is a gauge theory with couphng (3 = (3}. When N ^ oo 
at fixed this becomes a U{1) gauge theory. In D = 3 + 1 this theory has a strong 
couphng confining phase that is separated from a weak couphng Coulomb phase by a phase 
transition (probably first order) at [3c = 0(1) and a 'freeze-out' transition at P'^ = O^N"^) 
where fiuctuations between neighbouring elements of the centre become improbably small. 
These phase transitions will extend from = cxd to some finite f3a- In addition numerical 
calculations suggest that there is a first order line that crosses the /?/-axis (for iV > 5, at 
the bulk transition) and also the /3a-axis (apparently for all N) and to which the first order 
line from (3c is attached. This phase diagram has been explored in some detail for SU{2) J2| 
and SU{3) 18j One can interpret this phase structure in terms of the condensation of Zn 
monopoles and vortices These involve large plaquette values and so the first order bulk 
transition involves a large jump in the average plaquette. It appears, not surprisingly, that 
the would-be third-order Gross- Witten transition is subsumed in this jump |15j . 

By contrast, in the analytically tractable case of D = 1 + 1, the limAr^oo^Af ~ 
theory aX Pa = oo will have no finite-/^ phase transition by the same arguments used in jT] 
for SU{N < oo) gauge theories. One can also readily show that the finite-A^ cross-over along 
the f3f axis, which becomes the Gross- Witten transition at iV = oo, is matched by a similar 
cross-over and N = oo transition along the (3a axis. Indeed using the same change of variables 
and notation as in jl] one sees that Z{Pf = 0,Pa) = where V = LqLi and 



dWexpi^TifW^TifW} 



1 _ ftL 



(9) 



Here we have expanded the exponential and then used eqn(41) of [I] which is valid for —>■ oo. 
(This argument is casual with limits and is at most valid coming from the strong coupling 
side [T].) Now 

{TTaU, + l) = {\TTjU,\') = {\TlfW\') 

^ 1 



1 



Af2 



(10) 



Clearly as we increase from strong coupling there must be some non-analyticity in (Trat/p) 
at or before the value (3a/N'^ = 1- The corresponding non-analyticity in Z = represents a 
phase transition. For a precise derivation we refer to PD| . 

In D = 2 + 1 the limAf_^oo Zn ~ U{1) gauge theory at = cxo has no phase transition at 
finite (3 and is linearly confining at all (3 due to the the plasma of f/(l) monopole-instantons. 
In the Zat theory there is freeze-out transition at some (3 = 0{N'^) |21j. Numerical calculations 
in SU{2) suggest that the line of phase transitions descending from this point into the 
finite Pa plane ends before reaching the /3a = axis, so that there is only a peak in the specific 
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heat, but no phase transition, on the (3f axis. So all this suggests - albeit on limited evidence 
- that at finite N the D = 2 + 1 {f5a,f5f) phase diagram contains only smooth cross-overs 
(except for the freeze-out transitions), much like D = 1 + 1 and in contrast to the finite- 
phase transitions in = 3 + 1 . In this paper we shall show that this is indeed the case along 
the axis. 

3.4 Finite volume transitions 

Consider a. D = 3 + 1 SU (N) gauge theory on a LqL^ lattice where Lq ^ L. If we increase 
/3 from small values, then we will encounter a deconfining transition at a{f3)LQ = l/T^ (first 
order for > 3 and second order for = 2 [HI |3] ) . A convenient order parameter for this 
transition is the Polyakov loop, (/^=o), which acquires a non-zero expectation value in the 
deconfined phase. The transition is a crossover for finite L and sharpens to a phase transition 
as L ^ oo. As A^ t the transition becomes sharper on ever smaller volumes jHlEj. so that as 
N (yo one will have a deconfining phase transition for L = Lq + e where e > is as small 
as we like. By continuity one would expect an A^ = cxd phase transition even as e — > 0, i.e. on 
an lattice. (Finessing any subtleties about orders of hmits.) 

This makes contact with calculations [ZllHl that have shown that if we are on a lattice and 
increase (3, then there will be a crossover, sharpening to a phase transition at A^ = oo, which 
is characterised by one of the Polyakov loops, (/^), with fj, chosen at random, acquiring a non- 
zero expectation value. In fact one finds [HIHl that this is only the first of a sequence of phase 
transitions. At a second higher value of (3 there is a second phase transition where another 
Polyakov loop, (l^) with v ^ ^ again chosen at random, acquires a non zero expectation value. 
And there is some evidence that as (3 is increased further there are similar transitions along a 
third and fourth direction [71 |H]. Moreover one finds [71IH] that these transitions are associated 
with a gap forming in the eigenvalue spectrum of the appropriate Polyakov loop, just as one 
finds for the plaquette at the Gross- Witten transition. It is therefore natural to think of these 
finite-volume transitions as being in some sense strong-to-weak coupling transitions. 

As we argued above, the first of these N = oo transitions appears to be nothing but the 
N = oo deconfining transition and should occur at /? = Pco where 

a(/3,jLo = 1/Tf =^ (11) 

One can make analogous arguments for the existence of a sequence of transitions on an 

lattice with Lq <^ Li <^ L2 <^ L^. To see this, consider the following steps. 
(1) Increase (34^ = j3 to very high temperatures, T = l/a(/5)Lo ^ Tc. In this regime we will 
have the familiar dimensional reduction to an effective D = 2 + 1 SU {N) gauge theory coupled 
to adjoint scalars that are the time-translationally invariant remnant of the Aq gauge field 
|2Sj • To leading order the gauge coupling and mass of the scalar of the effective D = 2 + 1 
gauge-scalar theory are [21] 

9l = 9i{T)T ; ml<xgl{T)T'. (12) 

So ma/gl = 0(1/ g4{T)) and at high enough T the D = 3 + 1 gauge theory reduces to the 
SU{N) gauge theory in D = 2 -f 1 on a Li -C i^2,-^3 lattice. As we increase P4 = P we 
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simultaneously increase = 2N/ag^ ~ P^Lq (neglecting the difference between gl{a ^) and 
gl{T)). This D = 2 + 1 gauge theory will deconfine at 

a(/?,jLi = l/Tf =3 (13) 

at which point (/^=i) acquires a non-zero vacuum expectation value. We can estimate the 
corresponding critical value of /?(= P^) to be 

P,, ~ OMN^^. (14) 

To arrive at this estimate we use (Tc/y/a)D=3 ~ 0.9 and y/a/g^N ~ 0.198 (13j, together 
with eqn()12j). For finite N this will be a cross-over, but we expect (for the same reasons as in 
one higher dimension) that as — oo one will have a phase transition on any volume where 
L2, = Li + e, for any fixed e however small. In the limit we thus expect the transition on a 
symmetric LqL^ lattice with Li ^ Lq. If we now reduce Li/Lq then we see from eqn flT^ that 
/9ci will begin to decrease. However long before its value reaches Pc^ the value of T will have 
become small enough that we cannot neglect the adjoint scalar (and its self-interactions) and 
our estimate in eqn (|T^ ceases to be useful. Nonetheless it appears plausible that the second 
N = 00 transition that has been observed on lattices [3 |Hj is the continuation as Li Lq 
of this three-dimensional deconfinement transition. 

(2) As we increase P beyond Pc^ on our Lq -C Li -C L2 -C L3 lattice, T^"^ = 1/aLi will 
become ever larger, and eventually the system will undergo a further dimensional reduction 
to a D = 1 + 1 SU{N) gauge theory with adjoint scalars [221 . To leading order the gauge 
coupling and scalar mass of this effective D = 1 + 1 gauge-scalar theory are [23 

gl = glT^=' ; ml ^ gllog{aT^='). (15) 

While the pure gauge D = 1 + 1 theory has no propagating degrees of freedom and is too 
trivial to deconfine, the presence of the adjoint scalars renders it a non-trivial confining field 
theory which we would naively expect to deconfine at some =2 = l/a{pc2)L2. We estimate 
the corresponding critical value of the coupling /3(= P4) to be 

P^^r^OASr'N'-^ ■ r = ^ (16) 

using the value a ~ {0.8fglTN/3 extracted from [22] • Because we are in Z) = 1 + 1 the 
high-T phase cannot have a true non-zero expectation value for {1^=2), as will be discussed 
more explictly when we come to our numerical calculations below. Again we would expect 
that at A^ = 00 this deconfining transition will appear on lattices with finite L3 and even as 
L3 —>■ L2. It is again plausible that this N = 00 transition (although not our estimate in 
eqn ()l(j|l ) will survive as L3 —>■ L2 ^ Li ^ Lq, thus making the connection with the third 
transition observed on L"^ lattices |3 (Hj . 

(3) If we increase P further we come to consider a field theory with a finite Euclidean time 
extent given by aL^ living in an infinitesimal spatial volume a^LoLiL2. Such systems can in 
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principle have deconfining phase transitions [211 although whether this one does or not we do 
not attempt to make plausible by a simple argument. If it does exist then it would provide 
the final step in our cascade of N = oo phase transitions Pco Pci ^ Pc-i on our 

Lq <^ Li <^ L2 L-3 lattices. 

The above discussion has taken the D = 3 + 1 SU{N) gauge theory as its starting point. 
It is obvious that we could equally well have started with a Lq -C Li -C L2 15 = 2 + 1 SU{N) 
gauge theory and followed that through a cascade of deconfining N = 00 transitions. 

4 Results 

4.1 Preliminaries 
4.1.1 Phase transitions 

At a phase transition appropriate derivatives of the partition function Z diverge or are dis- 
continuous. (Strictly speaking of ^logZ where V is the volume.) The lowest order of such 
a singular derivative determines the order of the phase transition. For Z or its derivatives to 
be singular, we require an infinite number of degrees of freedom, and this usually demands 
an infinite volume, with a cross-over at finite V sharpening to the appropriate singularity as 
\^ 00. As ^ CX3 we have the possibility of a new kind of phase transition that takes 
place in a finite volume with the infinite number of degrees of freedom being provided by A^. 

With the standard plaquette action, a conventional first order transition has a discontinuity 
at V = 00 in the average plaquette, 

{up)=N~'d\ogZ/dp (17) 

where Np is the number of plaquettes. (We denote the space-time volume by V and note that 
V = Np/3 in D = 2 + 1.) At finite V this discontinuity is a rapid crossover so that 

d{up)/df3 = Np^d^ log Z/df3^ = C (18) 

diverges at the critical coupling (3 = Pc ^ Np —> 00. (Here C is the specific heat.) This 
divergence is linear in V since the cross-over between the two distinct values of the plaquette 
occurs in the range of f3 — f3c where there is back-and-forth tunnelling and this range is 0{1/V) 
as we see from the linear approximation 

AFi/3) = VAfifS) (X (/5 - /3,)V ~ 0(1) (19) 

to the free energy (density) difference, AF, /, between the two phases. 

A conventional second order transition has a smooth first derivative of Z but a diverging 
second derivative and specific heat C — 00 as V ^ 00. Defining to be the average value 
of Up over the space-time volume for a single lattice field, we easily see that the specific heat 
can be written correlation function: 

C = Np{{u^-{u;)f) = Np{{u;^)-{u^)^) 

= Y,{{up-{up)){up,-{up))) (20) 
p 
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where Pq is some arbitrary reference plaquette. It is clear from eqn(j2nj that the divergence of 
C a.s Np —>■ oo imphes that there is a diverging correlation length - the standard signal of a 
second order phase transition. 

A conventional third order transition has smooth first and second order derivatives but a 
singular third-order derivative, C = N~^d^ log Z/df]"^, at V = oo. This may be written as 

C = -^ = Nl{{u;-{u;)f) 

= N'piiu^') -3{u^){u^') +2{u^)'). (21) 

Note that if the fluctuations of were symmetric around (u^), as they are for (3 = 0, then C 
would be zero, so 

limC"(/5) = 0. (22) 

It should be clear that the higher the order of the transition, the greater is the statistics 
needed to determine its properties to a given precision. In particular, identifying third-order 
transitions is already a formidable numerical challenge, and we do not attempt to look for 
transitions that are of yet higher order. 

Since we are particularly interested in transitions that develop as — oo and since we 
know that, in general, fluctuations in the pure gauge theory decrease by powers of in the 
large- A^ limit 1121 HZ] it is convenient to define the rescaled quantities 

C2 = N^ xC ; C3 = N^xC' (23) 

which one expects generically to have finite non-zero limits when N —>■ oo. (Note that the 
increasing power of A^ simply matches the increasing power of d/djS = 2N'^d/d'j, where 7 is 
the inverse 't Hooft coupling defined in eqn(jH).) The signature of a phase transition which 
is only present for N = 00 will be a crossover for finite A^ at which fluctuations decrease 
more slowly than the naive power of l/N"^. If, therefore, we find a crossover in C2 or C3 
which does not sharpen with increasing volume at fixed A^, but rather becomes a divergence 
or a discontinuity only in the large- A^ limit, then this will indicate a second- or third-order 
A^ = 00 phase transition respectively. (Provided of course that {up) is continuous so that 
there is no large- A^ first order transition.) 

Large- A^ phase transitions can have an unconventional behaviour. Consider for example 
a second order transition characterised by a value of C2 that diverges at some A = Ac as 
A^ ^ 00. This may indeed be due to a correlation length ^ that diverges (in lattice units) 
as A^ ^ 00: .^(Ac) oc A^" ; a > 0. However there is another, less conventional, possibility: 
the correlation length may be finite and it may be that local plaquette fluctuations have an 
anomalous A^-dependence at the critical point: (Up) / (up)"^ — 1 oc A^"~^ ; a > 0. 

Since large- A^ phase transition can arise from fluctuations that are completely local - as 
in D = 1 + 1 where the lattice partition function factorises - it is also useful to consider local 
versions of the quantities C2 and C3: 

P2 = N^x ((V) - {upf) (24) 
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and 

P3 = iV^ X ((n/) - 3K)(V) + 2(«p)'). (25) 

These are the contributions to C2 and C3 from fluctuations of individual plaquettes, i.e. 
neglecting correlations between plaquettes. So calculations of P2 and P3 require much lower 
statistics than C2 and C3 to achieve the same accuracy, and this will be particularly useful 
at the largest values of A^. Of course, divergences or discontinuities in P2 or P3 will normally 
imply divergences or discontinuities in C2 or C3, even if the latter are not visible in the 
statistical noise of the numerical calculation. Note that in 1+1 dimensions, where the theory 
factorises and there are no propagating degrees of freedom, we have P2 = C2 and P3 = C3 at 
any A^. 

The eigenvalues of an SU{N) matrix such as the plaquette, are gauge-invariant, and so we 
can use them to gain additional information to that encoded in the above correlators of low 
powers of traces of plaquettes. Indeed, at the D = l-\-lN = oo Gross- Witten transition a 
gap opens in the eigenvalue density of the plaquette pp. That is to say, while in the strongly- 
coupled phase the eigenvalue density is non-zero for all angles — tt < a < tt, in the weakly- 
coupled phase it is only non-zero in the range —ac<<y<ac, where < vr. To search for a 
similar transition in 2+1 dimensions we will measure the total plaquette eigenvalue density. 
Using the eigenvalue density directly to search for a gap is difficult since for any finite N there 
is not a true gap but instead a (near-)exponential drop in the eigenvalue density as a — >• ivr, 
so very good statistics are required to observe changes in the exponentially suppressed tails of 
the density. To avoid relying solely on the accurate calculation of these tails, we also calculate 
the fluctuations of the eigenvalues around their average values, (A^) — (Aj)^, where Aj is the 
zth eigenvalue when ordered by its phase. (Recall that the eigenvalues of SU{N) matrices 
are pure phases Xj = exp{iaj}.) In particular, we shall calculate the normalised ratio for the 
extreme (smallest) eigenvalue: 



2 



(for N even) . This is motivated by the situation in the N 00 limit in D = 1 + 1 where the 
eigenvalue density at the Gross- Witten transition, 7 = 1/2, is P 2^(1 + cosa). This density 
approaches zero as a ^ ±tt but is finite at a = 0, so we expect the fiuctuations to be 0(1/A^) 
in Ajv while they can be up to 0(1) for Ai. Thus we expect Rp to diverge at the Gross- Witten 
transition, and it may provide a useful observable in our search for a similar transition in 
D = 2 + l. 



4.1.2 Wilson loop non-analyticities 

To investigate the possibility that Wilson loops undergo some analogous non-analyticity as 
their area passes through some critical value, Acru, we calculate Wilson loops of a fixed size, 
rii X n2, in lattice units and increase f3 so as to decrease the lattice spacing a and hence the 
area, A = ani x an2, in physical units. If there is a non-analyticity at A{f3c{'ni,n2)) we can 
then vary rii, n2 so as to check whether the transition occurs at a fixed area in the continuum 
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limit, when expressed in units of say g^N , i.e. whether 



hm f ) (27) 



is finite and non-zero. Since all the evidence is that the D = 2 + 1 SU(A^) lattice gauge 
theory has no phase transition, at zero temperature, once A is on the weak coupling side of 
the bulk transition, we expect any Wilson loop non-analyticity not to correspond to a phase 
transition of the whole system. This will be an important constraint on what are the important 
observables to calculate. We also expect that any such transitions will be cross-overs at finite 
A^, becoming real non-analyticities only at = oo. This is because we can imagine that they 
are driven by the degrees of freedom close to the critical length scale, and that we need these 
to be infinite in number for a real non-analyticity. 

We remarked in Section IT^ that a non-analyticity in the eigenvalue spectrum of the Wilson 
loop is known to occur |9, in the D = l + lN = oo gauge theory at the critical area given 
in eqn((7)). We have performed numerical lattice calculations in this theory for large and 
find that the lattice critical area is very close to the continuum one for Wilson loops that are 
2 X 2 or larger. To be more precise let us denote the product of link matrices around a square 
nx n Wilson loop by [/^^" and its trace by u'^^^ = -^ReTrlf/^^*^} which we generically write 
as Uw Then we find that the non-analyticity occurs when reaches a particular value 

{u^) ^ e-\ (28) 

As an example we show in Fig. ^ the eigenvalue spectrum of a 3 x 3 Wilson loop in D = 1 + 1 
for A^ = 48 at A = 0.7971 where the trace satisfies eqn ()28|l and we compare it to the continuum 
expression obtained from [01 Elj We clearly have a very good match (apart from the A^ = 48 
bumps that arise from the eigenvalue repulsion in the Haar measure). Now we know that in 
D = 1 + 1 the Wilson loop factorises into a product of plaquettes 

= {up)^ (29) 
and that {up) = 1 — A/4atA^ = oo 1^. Putting all this together, we have 

1 - - j = 4^ ~ e 3^ ~ e (30) 

which we observe is nothing but the continuum relation in eqn((Zj). These numerical calcula- 
tions show that lattice corrections are small except for loops smaller than 2x2, such as the 
plaquette that has its non-analyticity at the Gross- Witten transition where A = 2. Because 
of the factorisation in eqn()29p the trace of will be analytic in the (bare) coupling when 
this gap in the eigenvalue spectrum forms (except for the very smallest loops where it occurs 
at the Gross- Witten transition) and so it is not immediately obvious what is the significance 
of this gap formation. What this tells us, nonetheless, is that we should not only search in 
D = 2 + 1 for non-analyticities of traces of Wilson loops, but also for such eigenvalue gap 
formation. 
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Our initial question will be whether (u^) has a non-analyticity in A for some value of the 
area, A/a"^. To investigate this we can look at (u^) and its derivatives as a function of the 
bare coupling A. The derivatives can be expressed as correlation functions in the usual way 
e.g. 

d{uj)/d(3 = Np{{uwU^) - {u;){uyj)) 

= ^{{U^- {Uw)){Up~ {Up})). (31) 

p 

Because the whole system has no phase transition we expect that the non-analyticity will be 
visible in the 'local' correlators analogous to those in eqns ()24f25|l e.g. 

= N^x {{u^up) - {u^){up)) (32) 

where Up is the average of the plaquettes that tile the minimal Wilson loop surface. We define 
in analogy to P3 in the same way. 

It is also possible that some non-analyticity might be present in just the fluctuations of 
Wilson loops rather than in the derivatives with respect to the coupling. Thus we also consider 
the correlators 

^nxn = ^ («^"') - «^")') (33) 

and 

^nxn = ^ (2«>^")3 - 3«^")(m:^^"') + (m^^"')). (34) 

that represent an alternative generalisation to Wilson loops of the quantities P2 and P3 defined 
in eans (l24l25jl . 

In searching for a gap formation in the eigenvalue spectrum of an n x n Wilson loop, we 
define the quantity i?"^" in direct analogy to the quantity Rp defined in eqn fj26|) . However 
this quantity is only useful for the bulk transition because the eigenvalue spectrum approaches 
the gap with a finite slope. We know that for larger Wilson loops in D = 1 -(- 1 [01 El the 
approach is with a diverging slope (at N = 00), and that R^^^ is not a useful observable 
in that case. We shall in fact find it much more useful to match the eigenvalue spectra in 
D = 2 + 1 and D = 1 + 1. That this is in fact possible is one of our most interesting results. 

4.2 Bulk transition 

In 3+1 dimensions the bulk transition is easily visible as a large discontinuity in the action 
for > 5 (where the transition is first order) and as a (finite) peak in the specific heat for 
N < 4 (where the transition is a crossover). We have searched for an analogous jump or rapid 
crossover in 2+1 dimensions, in particular around 7 = j3/2N'^ ~ 1/2. 

In FigI21we display the values of the average plaquette for SU{Q), SU {12), S'f/(24) and 
SU{4&) as obtained on 6^ lattices. At 7 ~ 1/2 an L = 6 lattice has a size La^o ~ 3 and 
so is large enough that it should display a very sharp cross-over for a conventional first-order 
transition. This should be more so as t and (most) finite volume effects disappear. As a 
check we have repeated our calculations on 12^ lattices for S'f/(6) and have found no volume 
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dependence. What we see in Fig. El is that the action appears to be approaching a smooth 
crossover in the large- hmit, with no evidence for a first order phase transition either at 
finite or at = oo. 

Our results for the specific heat C2 for SU (6) and SU {12) are shown in Fig.lHl (For SU (24) 
and 5'[/(48) our accuracy is insufficient to get useful results for C2.) There is a clear peak 
around 7 ~ 0.42 which appears to be growing stronger with increasing N . For SU{&) we 
repeated our calculations on 12'^ lattices and found no volume dependence. This tells us that 
we are not seeing a conventional second-order phase transition at fixed N for which the specific 
heat peak grows as the volume increases (since a larger volume can better accommodate the 
diverging correlation length). So if there is a second-order phase transition here it would 
appear to be not at finite A^, but only at A^ = cxd. 

To search for a possible third order transition we calculate C3, but our calculations are 
not accurate enough to produce anything significant, even for S'f/(6). 

To improve our reach in A^ we calculate the quantities P2 and P3 defined in eqns (j24l25p . 
These represent the contributions to C2 and C3 made by the fluctuations of individual pla- 
quettes and are the quantities that reveal the Gross- Witten transition in D = 1 + 1. In Fig. |3] 
we show the values of P2 obtained for S'f/(6), SU {12), SU{24^) and SU{4S>). We observe, as 
expected, a dramatic reduction in the statistical errors as compared to C2 in Fig. 13 enabling 
us to look for fine structure in the /5-dependence. There is no significant evidence for a peak 
in P2 which indicates that if there is a second order transition at A^ = cxd, as suggested by 
the peak in C2, it will primarily involve correlations between different plaquettes rather than 
arising from the fluctuations of individual plaquettes. What we do see in P2 however is defi- 
nite evidence for a cusp developing with increasing A^, at 7 ~ 0.43, where the derivative of P2 
will suffer a discontinuity. This corresponds to a third-order transition at A^ = 00, just like 
the D = 1 + 1 Gross- Witten transition pQ. It is therefore useful to compare the D = 1 + 1 
and D = 2 + 1 cases in more detail. For this purpose we show in Fig. some numerically 
calculated values of P2 in -D = 1 -|- 1 SU{Q), SU{12) and SU{2A) gauge theories. (Recall that 
in D = 1 + 1 the factorisation of the partition function implies that C2 = P2-) We also plot 
the analytic results for SU {00) [1/. 



Apart from a small relative shift in 7 the results for D = 2 + 1 and D = 1 + 1 are remarkably 
similar, strengthening the evidence for a third-order N = 00 transition. 

To investigate this further, we show in Fig. IHlour results for P3 for SU{6), SU{12), S'f/(24) 
and SU{4:8). There is clearly an increasingly sharp transition as A^ increases around 7 ~ 0.43. 
For comparison we show in Fig. [7| corresponding numerical results for D = 1 + 1 (where 
C3 = P3) together with the analytic result for SU {00) 0: 



which has a discontinuity at the Gross-Witten transition at 7 = 1/2. It is clear that once 
again the the behaviour in D = 2 + 1 is remarkably similar to that in 1-1-1 dimensions. 




(35) 




(36) 



13 



We see further evidence for a Gross- Witten-like transition in our results for the ratio Rp 
defined in eq. 1211 Our results for SU{Q), 5f/(12), SU{2A) and SU{A8) are plotted in Fig. E 
There is a clear peak around 7 ~ 0.43 whose height increases rapidly with N, indicating that 
the fiuctations of the extreme eigenvalues are becoming much larger than fluctuations of the 
'middle' eigenvalue (the one nearest a = 0). Very similar behaviour occurs in 1+1 dimensions, 
for which our results are shown in Fig. El Indeed, in the iV — > 00 limit in 1+1 dimensions we 
expect that Rp will diverge, as discussed below eq. EHl This appears to be exactly what we 
see in Fig. (HI 

Finally we compare the D = 2 + 1 and D = 1 + 1 transitions directly by comparing 
the eigenvalue densities across the transition. We do this for SU(12) in Fig. El The eigen- 
value densities both below and above the transition are clearly very similar in 1+1 and 2+1 
dimensions. 

All the above confirms that D = 2 + 1 SU(A^) gauge theories possess an = 00 third- 
order strong-to-weak coupling transition that is remarkably similar, even in its details, to the 
D = 1 + 1 Gross- Witten transition. 

Despite this striking similarity, when we look in more detail we also observe significant 
differences between the bulk transition in 2+1 dimensions and the Gross-Witten transition. 
Comparing Figs. El and we see that there is a peak in the specific heat in D = 2 + 1 
which is simply not present in D = 1 + 1. From Fig. jH it is clear that this peak does not 
come from fluctuations of individual plaquettes, but must come from correlations between 
different plaquettes. To investigate this we consider the following particular contributions to 
the specific heat C2: the contribution from correlations between a plaquette and its neighbours 
in the same plane, which we label Cj; the contribution from correlations between a plaquette 
and its neighbours which share an edge but are not in the same plane, Co] and finally Cf, 
the contribution from correlations between a plaquette and the plaquettes facing it across an 
elementary cube. We include a factor A^^, as for C2. We find a clear peak, growing with A^, 
in our results for Co, plotted in Fig. [TTl The peak accounts for about half of the difference 
between C2 and P2. There is also a much weaker peak in Ci, approximately a factor of 15 times 
smaller, which also clearly grows with A^, at least up to A^ = 24. (The weakness of the signal 
means that we lose statistical significance for larger A^.) For Cf, where we happen to have 
results only for SU(6) and SU(12), we see in both cases a clear peak. This is almost exactly a 
factor of four lower than the corresponding peak for Co- Since each plaquette has four times as 
many out-of-plane neighbours as it has neighbours facing it across an elementary cube, this 
shows the correlation of a plaquette with its individual out-of-plane neighbours is in fact the 
same as with a facing plaquette. By contrast the correlation with the 'nearer' neighbouring 
plaquettes that are in the same plane (as measured by Ci) is very much weaker. This pattern 
is precisely what one would expect if the correlations were due to a flux emerging from the 
cube symmetrically through every face, i.e. due to the presence of monopole-instantons. 

If such monopoles are present, we would expect the correlation of the plaquette with itself, 
P2, to be also affected. These correlations of the plaquette with itself should be as large as 
with each of its eight out-of-plane neighbours, so this contribution to P2 should about one 
eighth of Co- Even for SU{4:8) this is only ~ 0.04, which is easily consistent with our results in 
Fig.m Of course, if there really is a second-order phase transition at A^ = 00, then eventually 
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we would expect to see a growing peak in P2. It is of course possible that there is no second- 
order phase transition at = cxd, but only a rapid cross-over, so that Co asymptotes to a 
finite value, and in that case there would not need be a pronounced peak in P2. This however 
seems a rather articial scenario. 

Since a second-order transition is usually associated with a diverging correlation length, 
we also measured the mass of the lightest particle that couples to the plaquette, in both SU{6) 
and SU (12). We calculate an effective mass at a separation of n lattice units: 



where 0o is the trace of a plaquette and 0„ is the trace of a facing plaquette lying in the same 
plane n lattice spacings away. This is not a zero-momentum correlator, so it will overestimate 
the masses. (We do not have a statistically significant signal from zero-momentum correlators.) 
Our results for ameff(l) are plotted in Fig. [T21 (Our results for n > 2 do not have a useful 
statistical accuracy.) We observe a dip in the effective masses near the transition, which 
becomes more significant for SU (12), particularly when we take into account the expected 
weak-coupling scaling behaviour, am oc 1/(3. While this is certainly consistent with a second- 
order cross-over, the masses are large, and if the correlation length is going to show any sign 
of diverging it is clear that it will be at a much larger value of than are accessible to our 
calculations. 

4.3 Wilson loops 

4.3.1 Traces and correlators 

In Fig El we show how (uy,) varies with A in some sample calculations. We see no sign of 
any singularity developing in this quantity, or in our simulataneous calculations of d{uw)/d\, 
in contrast to the growing peak we saw for C2 oc d{uw)/dX in Fig. El The local version of 
the correlator corresponding to the first derivative defined in eqn fl32j) . P^, is more accurately 
calculated and its variation with A is shown in Fig El and shows no evidence of a developing 
cusp that would suggest an = 00 singularity in the second derivative. Thus, at this level 
of accuracy, we see no evidence for any N = 00 non-analyticity in the variation of (uw) as a 
function of the coupling A. 

Given our uncertainty in the type of analyticity that might occur we have also considered 
it worthwhile to look at the quantities IV^^^" and IV^^" which are defined in eqs. OHl and 
El and which are alternative analogues of the quantities P2 and P3 for the plaquette. Some 
results for SU{Q), 5f/(12), SU{2A) and SU{A8) are plotted in Fig. d While there is a range 
of 7 = 1/A over which W2^'^ stops being constant and starts to decline, this 'transition' does 
not become sharper with A^, in contrast to the behaviour of the plaquette equivalent, P2, in 
Fig. m We see the same behaviour for 3x3 and 4x4 Wilson loops, only shifted to higher 
7, with no sign of the transitions becoming sharper as A^ increases. Our results for W^^"^ for 
SU(6) and SU(12) are plotted in Fig. El While there is a transition region in which W^^"^ 
becomes negative, just as one sees for P3 in Fig. (HI the transition does not become sharper as 




(37) 
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increases, unlike P3. Again we see the same behaviour for W^^^ and W^^^, with no sign 
of the transitions becoming sharper as N increases. 

All these results are in fact essentially identical to those we obtain in similar calculations in 
D = 1 + 1, where we know that the (uu,) is analytic in A except at the Gross- Witten transition. 

Finally we recall that for the plaquette the Gross- Witten transition is characterised by 
a divergence in the relative fluctuation of extremal eigenvalues, as shown in Fig. |H1 For 
n X n Wilson loops the analogous quantity, K^^^, shows no such behaviour, as we see, for 
the example of i?2x2 for SU{6), SU (12) and SU{2A) plotted in Fig. El This is perhaps no 
surprise, given that in D = 1 + 1 the eigenvalue gap formation for Wilson loops larger than 
the plaquette does not involve growing fluctuations of the extremal eigenvalues. 

4.3.2 Matching eigenvalue spectra 

Although we have found no evidence that the trace of a Wilson loop is non-analytic in A at 
some critical area, it is possible that there are more subtle non-analyticities of the kind that 
exist in D = 1 + 1 and which are associated with a gap forming in the eigenvalue spectrum. 

To determine numerically whether at some given A the spectrum p{a) in some region close 
to a = ±7r will extrapolate exactly to zero when A^ ^ 00 is clearly a delicate matter, given 
that the values at finite A^ from which we extrapolate are already extremely small. 

So to search for such non-analytic behaviour we explore the strategy of directly comparing 
Wilson loop eigenvalue spectra in 1+1 and 2+1 dimensions. We first evaluate the spectrum 
in 1+1 dimensions at the critical coupling at which the gap forms. A true gap only forms at 
N = 00; for finite A^ we use the same value of the critical 't Hooft coupling, ITI]] 

Ac = - = 4(l-e"x^), (38) 

7c 

where A is the area of the Wilson loop in physical units. At this coupling the expectation 
value of the trace of the Wilson loop is, using eqn® P ITU]. 

(u^) = {{u,)}^^^' i^l - =e-\ (39) 

which is the same value as at the critical coupling in the continuum limit. Note also that as 
a ^ and A/o? 00, eqn (jnH|) reduces to eqn((7j) as it should. Having obtained the spectrum 
(numerically) in = 1 + 1 for a given size Wilson loop (in lattice units) and for a given value 
of A^, we then calculate the eigenvalue spectrum in D = 2 + 1 for the same size loop and for 
the same A^, varying the coupling to a value where the two eigenvalue spectra match. 

We find that it is always possible to achieve such a match, for any A^ and for any size of 
Wilson loop. We show an example in Fig. El where we compare the eigenvalue density of the 
3x3 Wilson loop in SU(12) in 1+1 dimensions to the density in 2+1 dimensions, at a coupling 
chosen to give the best match. In Fig. Elthe coupling m D = 1 + 1 is Ac, the coupling at which 
the gap forms. The spectra are clearly very similar and indeed indistinguishable on this plot. 
We also find that the spectra can be matched when they are away from the critical coupling: 
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we show this in Fig. E[ where we choose a higher value of 7 to illustrate that the matching 
continues to work after the gap forms. In Figs. Eland ^Iwe also plot the analytically known 
spectra [01 . in the continuum limit of the N = 00 theory in D = 1 + 1 at the corresponding 
couplings. These clearly match the corresponding finite-A^ spectra very well, except in two 
respects: the latter have bumps which arise from the eigenvalue repulsion that is a well- 
known characteristic of the Haar measure, and the finite- spectrum is not precisely zero in 
the region of the 'gap'. 

The fact that at finite but large we can match so precisely the D = 1 + 1 and D = 2 + 1 
eigenvalue spectra for couplings at and above the D = 1 + 1 transition, provides convincing 
evidence that the Wilson loops in the D = 2 + lN = oo theory also undergo a transition 
involving the formation of a gap in the eigenvalue spectrum. 

In Fig. 1201 we plot the eigenvalue spectra of 2 x 2, 3 x 3, and 4x4 loops in SU(6) in 
2+1 dimensions. The three couplings have been chosen so as to give the best match to the 
eigenvalue spectra of Wilson loops of the same size in D = 1 + 1 at Ac- We see that the 
three spectra are essentially identical. Moreover the critical value of the D = 2 + 1 coupling 
7c = 1/Ac appears to grow linearly with the size of the L x L loop, suggesting that there is a 
finite critical area for gap formation in the contimuum limit: ITU] 

XlL' {ag'NYL' = (/Ar)^^,^,. (40) 

As we shall see below, in Section [4. 3. 4| this is nearly but not quite the case. 

It turns out that all the above is an immediate corollary of a much stronger and rather 
surprising result concerning the matching of Wilson loop eigenvalue spectra in 1+1 and 2+1 
(and indeed 3+1) dimensions. 

The general statement is that if take an nxn Wilson loop ?7^^" in the SU(A^) gauge theory 
and calculate the eigenvalue spectra in D and D' dimensions, we find that the eigenvalue 
spectra match at the couplings and A^' at which the averages of the traces m^^" = 
^ReTrjt/;;;^"} are equal: 

(«r"(AD))D = {uT-i\D'))D'. (41) 

We have tested this matching for D = 1 + 1 and D = 2 + 1 over groups in the range N = 2 
to A^ = 48 and for Wilson loops ranging in size from 1x1 (the plaquette) to 8 x 8 and, in 
2+1 dimensional, for couplings from A = 4.0 to A = 0.40. We have in addition tested it in the 
deconfined as well as in the confined phase. Some sample calculations in D = 3 + 1 have also 
been performed strongly suggesting that the same is true there. 

The fact that such a precise matching is possible implies that the eigenvalue spectrum is 
completely determined by A^, the size of the loop, and its trace. Hence the eigenvalues are not 
really independent degrees of freedom, which is unexpected. Moreover we have seen in Fig. EOl 
a demonstration of the fact that the spectra of Wilson loops that are 2 and larger can also be 
matched with each other. The matching occurs at values of the traces that are the same as 
those in D = 1 + 1 where they are calculable. In this sense, the size of the Wilson loop is not 
really an extra variable here. Finally, the A^ dependence is weak, and consists mainly of the 
two differences noted earlier. 
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Finally we remark that our results at this stage rely on a comparison that is visual and 
impressionistic. Ideally one would like to match the spectra by varying A continuously and 
this can be done, from nearby calculated values of the coupling, by standard reweighting 
techniques. In addition it would be useful to quantify any differences (which must be very 
small) with a standard error analysis. We intend to provide such analyses elsewhere 

4.3.3 Polyakov loops 

We have also investigated the eigenvalue spectra of Polyakov loops as defined in Section El 
These are products of link matrices that wrap around one of the space-time tori (and are 
of minimal length unless specified otherwise) i.e. they can be thought of as non-contractible 
Wilson loops. They provide the conventional order parameter for the deconfinement phase 
transition. As one crosses this transition the Polyakov loop that winds around the time (tem- 
perature) torus acquires a non-zero expectation value. This corresponds to the spontaneous 
breaking of a global centre symmetry in the Euclidean system. To simulate the system at 
temperature T we use a L'^Lq lattice with Lg ^ Lt so that T = I/oLq. As N grows one can 
weaken the inequality, so that one can take Lg ^ Lf as N ^ oo while still maintaining the 
thermodynamic interpretation and the sharp phase transition. 

We calculated the eigenvalue spectra of timelike Polyakov loops in SU (12) on L'^Lq lattices. 
We found that it is always possible to match the Polyakov loop eigenvalue spectra to those 
of Wilson loops in 1+1 dimensions (and hence also to Wilson loops in 2+1 dimensions) by 
choosing couplings at which the trace of the Polyakov loop equals that of the Wilson loop 

\{l^=o)\ = M (42) 

where, as we have seen, the size of the Wilson loop does not matter to a very good approxi- 
mation. (We take the modulus because the Polyakov loop is proportional to some element of 
the centre in the deconfined phase and the modulus effectively rotates that element to unity. 
The eigenvalue spectrum also needs to be rotated by the same centre element.) This matching 
has the corollary that the Polyakov loop eigenvalue spectrum will develop a gap at = oo 
when its trace crosses the critical value |(/^=o)| = For A^ > 4 the deconfining transition at 
T = Tc is strongly first order and the value of | (/^=o) | will jump from | {lfj,=o) | = at T < Tc to 
some non-zero value for T = T+. The latter value will typically be greater than for small 
Lq, i.e. for coarse lattice spacings, and will — as a — *• and hence Lq — >■ oo. Moreover for 
fixed Lf) the trace increases with increasing T. (See Section [4.3.41 for why this is so.) Thus for 
coarse lattice spacings we expect the gap formation to occur at the phase transition, T = Tc, 
while for larger Lq it will not coincide with the deconfining transition; instead it will occur at 
some T > Tc. The critical value turns out to be Lq = 7. Thus in the continuum limit the gap 
formation in timelike Polyakov loops does not occur at T = but rather at T = oo. 

As a numerical example of the eigenvalue matching we show in Fig |22l the eigenvalue 
spectra of the timelike Polyakov loop just below and just above the deconfinement transition 
for Lt = 4, together with a 3 x 3 Wilson loop spectrum in 1+1 dimensions at a coupling chosen 
to match the spectrum of the deconfined Polyakov loop. The spectra clearly match closely. 
Since the 1+1 dimensional 7 is above 7c for the 3x3 loop, the Wilson loop will develop a gap 
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at this coupling in the large- limit. Hence the Polyakov loop will presumably also develop 
a gap. 

Finally we recall that as — oo the deconfining transition occurs on smaller spatial 
volumes Lg Lq so that at = oo one can discuss the transition on a lattice. Taking 
into account the fact that our preliminary results |15] indicate that all the above carries over 
to Wilson and Polyakov loops in D = 3 + 1, we can make direct contact with the observation 
in [71 IH] that on an lattice the Polyakov loop develops a gap when it develops a non-zero 
expectation value. 

4.3.4 Theoretical interpretation 

The fact that at = oo there is a gap at weak coupling in the eigenvalue spectra of Wilson 
loops, has a simple explanation in the theory of Random Matrices. (See e.g. j2HI for a recent 
review.) At A^ = oo the Gaussian Unitary Ensemble (GUE) of complex Hermitian N x N 
matrices generates an eigenvalue spectrum that is the well-known Wigner semicircle 



In weak coupling, when /3 — > oo, the SU(N) link matrices can be expanded in terms of the 
Hermitian gauge potentials and it is very plausible that the averages involved in the calculation 
of Wilson loops fall into the same 'universality class' as the GUE. That is to say, once the 
eigenvalues of Wilson loops are clustered close to unity, the fact that the phases are on a circle 
rather than on the line becomes irrelevant and the phases (suitably rescaled by the coupling) 
should be distributed according to the semi-circle in eqn ()43j) . In fact this is precisely what we 
find. Thus the existence of a gap in the eigenvalue spectrum at weak coupling has a rather 
general origin in terms of Random Matrix Theory. 

On the other hand we know that in a confining theory 



which requires a nearly flat eigenvalue spectrum in [— vr, +7r]. Thus as we decrease the lattice 
spacing, the eigenvalue spectrum of a L x L Wilson loop must change from being nearly 
uniform to eventually having a Wigner semicircle gap. Thus at some bare coupling it must 
pass through a transition where the gap forms. 

For this gap to be physically significant, it must occur at a fixed physical area in the 
continuum limit. However, as we shall now see, this is not the case for either 2 + 1 or 3 + 1 
dimensions (in contrast to D = 1 -|- 1). The reason is the perturbative self-energy of the 
sources whose propagators are the straight-line sections of the Wilson loop. (Often referred 
to as the 'perimeter term'.) The leading correction is given by the Coulomb potential Vc{r) 
at the 'cutoff' r = a. For a Wilson loop whose size is / x / = aL x aL in physical units, this 
correction is 
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using the fact that Vc(r) oc g'^N logr, g'^N/r and A = ag^N,g'^NmD = 2+1, 3+1 respectively. 
Let us, for illustrative purposes, assume that the full potential is given by this self-energy and 
the area piece, aA = a^crL^, that comes from linear confinement. Then we have 

(m^) ocexpjcALlogA-c'A^L^} : D = 2 + l (46) 
using the fact that a'^a oc {ag'^NY = A^ and log a = log A + ■ ■ ■ in D = 2 + 1, and 

(m^) ocexpjcAL-c'e-'^L^I : D = 3 + 1 (47) 

using the fact that o?a oc exp{— c^/fi'^A^} in D = 3 + 1, where is given by the coefficients of 
the 2-loop renormalisation group equation. 

Consider first the D = 2 + 1 case in eqn ()4(jj) . Since XL = ag^NL = g'^Nl is the length 
scale in physical units, we see that if it were not for the weakly varying log A term in eqn ()46p . 
the Wilson loop trace would be the same on the lattice and in the continuum (up to the usual 
0{a?) lattice corrections). That is to say, we expect that as we approach the continuum limit, 
A ^ 0, the critical area for gap formation will vanish 

Acr^t OC — ^ ^ (48) 

(log A) 2 

rather than tending to some finite limit. At coarse a the logarithmic correction will be weak 
and one might well be tempted to perform an extrapolation to the continuum limit that does 
not include it. We illustrate all this with a numerical calculation of the coupling, and hence 
lattice spacing, at which Lx L loops develop a gap. We define the appearance of a 'gap' in our 
finite-A^ calculations as the coupling at which the spectrum is closest to the spectrum of the 
Lx L loop in 1+1 dimensions at the coupling 7c, for the same A^. For SU(2) we calculated this 
coupling for L up to 8 on 16^ lattices. For SU(6) we calculated up to L = 4 on 6^ lattices. We 
show our results in Fig.|^ together with a best fit to the SU(2) data which has the asymptotic 
behaviour in eqn()48|). The numerical data shows deviations from linearity which could either 
be interpreted as low L corrections to an asymptotic scaling behaviour 7 = A~^ oc L, or as a 
logarithmic violation of this asymptotic scaling. From our above analysis we know the latter 
to be the correct interpretation. 

In contrast to the anomalous behaviour we see when taking the continuum limit of Xc{A), 
the large- A^ limit is achieved rapidly and smoothly. To illustrate this we list in Table ^ the 
coupling for which the 3x3 loop develops a gap for A^ G [2,48]. The critical coupling is 
essentially constant from SU{Q) onwards, showing that are in the large- A^ limit. Indeed, even 
for SU(2) the corrections are small. 

In the case of D = 3 + 1 the self-energy diverges linearly and will normally dominate 
the trace for all A in the weak coupling region. Thus we expect Acrit oc up to logarithmic 
corrections from the running coupling, so that the gap formation occurs in the deep ultraviolet 
as we approach the continuum limit. 

In contrast, in D = 1 + 1 where the Coulomb potential is linear V^(r) oc g'^Nr, the 
self-energy term contributes at most a mere lattice spacing correction that vanishes in the 
continuum limit. 
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From the above discussion we see that the anomalous behaviour of Acru as a — > arises 
from divergent self-energy contributions. If the source had a finite mass, so that the propagator 
was smeared over some range (5r ~ we would evaluate the Coulomb self-interaction at 
r = rather than at r = a and hence would replace ALlogA — » ALlog/x in eqn ()4fjj) . and 
AL — > A//i in eqn ljTfjl . Assuming the univerahty of the gap formation persists for such loops, 
we would then expect them to form a gap at a value of A^ru that is finite in the continuum 
limit if we have chosen /i to be finite in physical units. The value of Acru will of course depend 
on the value of fi. 

Similar considerations apply to Polyakov loops. 

While the above considerations make plausible the universality aspect of the gap formation 
in Wilson and Polyakov loops, they do not explain our most striking result which is that the 
complete eigenvalue spectra can be matched across space-time dimension and loop size by 
merely matching traces. 

Finally, whether the gap formation, and the associated non-analyticity, has any significant 
physical implications is unclear. For that to be so one would require that the gap should form 
at a fixed physical area Acru in the continuum limit. As we have seen that is not the case in 
D = 2 + lorinZ) = 3 + l and is only the case in D = 1 -|- 1, where there are no propagating 
degrees of freedom and so no 'physics' in the usual sense. One can imagine regularising the 
divergent self-energies so that Acru is finite and non-zero in the continuum limit, but then it 
would appear to depend on the regularisation mass scale /i used. 

4.4 Finite volume 

In Section 13.41 we argued that = 3 + 1 SU{N) gauge theories on Lq ^ <^ -^2 ^ L-s 
lattices undergo a series of = oo phase transitions at Pco ^ Pci -C /3c2 ^ Pc^- These phase 
transitions are essentially deconfining transitions, a[(3c^)Li = 1/T^^^~^, in a dimensionally 
reduced theory. We argued that continuity, and vanishing finite size corrections at large A^, 
link these transitions to the N = oo phase transitions on L'^ lattices discussed in |21 El HI • 

The same argument clearly holds for SU{N) gauge theories on Lq <^ <C L2 lattices 
in D = 2 + 1 . Here we provide some (very) exploratory numerical results in support of this 
scenario. We do so on lattices with a less than asymptotic ordering, Lq < Li < L2. 

4.4.1 First transition 

The first transition is the usual deconfining phase transition when Li, L2 ^ 00. It is second 
order for SU(2) and SU(3), either second or first order for SU(4), and first order for N > 5 
|24j . Because the latent heat for A^ > 5 is oc A^^, the cross-over on a finite Lq < Li, L2 lattice 
will become a first-order phase transition at A^ = 00. All this is well-established and does not 
require further numerical confirmation in this paper. 
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4.4.2 Second transition 



To search for the second transition we simulate SU{12) gauge fields on a L0L1L2 = 2 x 4 x 40 
lattice over a large range of 7 = P/2N^. We calculate the Polyakov loop around the /x = 1 
torus, average it over the given lattice field, and take the modulus: This provides the 
conventional order parameter for a deconfining transition with the Li = 4 torus providing the 
(inverse) temeprature. We plot results for the average of this, at each value of 7 in 

Fig. 1221 We also plot values of the plaquette difference, ((mqi —^02)), which should also reflect 
such a transition. We see in Fig. |2S1 a very clear signal for a transition at 7 ~ 3.2 in both 
quantities. This occurs at a temperature Tc^ = T^^^^^ = l/4a(7 ~ 3.2) in the reduced theory. 
In units of the usual deconfining temperature of the D = 2 + 1 gauge theory, = T^^"^"*"^, 
this amounts to ~ 3Tco. 

The rapid, steep crossover suggests that the transition is first order. We find support for 
this when we plot a histogram of the values of |/i| at some of the values of 7 in the cross- 
over region. This is illustrated in Fig. I^for an ensemble of SU{12) fields on a 2 x 4 x 80 
lattice at 7 = 3.368. We see a clear two state signal, with the peak at small |/i| being 
naturally interpreted as coming from fields in the confining phase and the peak at large as 
coming from fields in the deconfined phase. Such a two-state signal is typical of a first order 
deconfining transition. 

The next question is whether this cross-over will sharpen into an actual phase transition 
in the two interesting limits: when we increase the spatial volume (here just aL2) at fixed 
N] or when we increase N at fixed volume. In addressing the former question we need to 
remark upon some special features of first-order transitions in the effective D = 1 + 1 theory 
that we are discussing here. The high temperature deconfined phase is normally characterised 
by a centre symmetry breaking so that / ~ c{(3) exp{i27m/N} where c{f3) is a self-energy 
renormalisation factor. Two such phases, characterised by n and n' say, can coexist and will 
be separated by a domain wall whose tension we expect [22] to be 



rp2 ^ 

Ok oc kiN - k) , ; A; = In - n'l , T = — . (49) 

In one spatial dimension the domain wall is just a 'point' and so the usual energy /entropy 
arguments tell us that at T = 1/ aLx > T^^ the field will break up into domains of typical size 

Arocexp| + ^} (50) 

Thus at any T if we take L2 ^ 00 the volume will consist of a 'gas' of domain 'walls', and 
hence domains, and on the average these will be equally distributed amongst all the centre 
phases, so that (|r^=i|) — >• 0. However on volumes that satisfy Li -C L2 -C Ar we will typically 
be in one domain and will thus have the usual deconfining signal of a non-zero value for |/^=i|. 
In addition it is clear that the lightest mass, rup, coupling to the = 1 Polyakov loop will 
not vanish at T > Tc^ but will approximately satisfy rrip oc exp{—cNT/y/X}. Note that this 
mass decreases with increasing Li = 1/aTi in contrast to the stringy behaviour, nip oc Li, in 
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the confining phase. It is clear from the above that our conventional signals for being in a 
deconfined phase become more complicated to interpret in D = 1 + 1. 

Similar considerations apply to the deconfining transition itself. Suppose that there are 
confining and deconfining phases that differ by a free energy density /. At T = we 
have / = so that the typical field will consist of a 'gas' of domain walls of typical size 
Ar oc exp{+(Tcrf/T} where ad is the free energy of the confining- deconfining interface. This 
is the essential difference with higher dimensions. For large enough volume (= L2) half the 
domains will be confining and half will be deconfining. Let us now increase the temperature 
T a little above Tc^. Then / = eo(T — TcJ near Tc^, where eo = [m]^ in D = 1 + 1. If T — Tcj is 
small enough, then Arf{T)/T <C 1 and the fraction of the volume that is still in the confined 
phase will be oc exp{— Ar/(T)/T} ~ 0(1). That is to say, the transition will take place over 
a range of temperatures AT that is no smaller than 

AT 

— oceo exp{-acd/Tc,} (51) 

-l-ci 

and this remains non-zero in the infinite volume limit. This implies that in D = 1 + 1 there 
cannot be an infinitely sharp first order transition even in the large volume limit. However, 
because both ad (probably) and / (certainly) grow oc iV^, there can be a phase transition at 
N = 00, and this can occur at finite volume. 

Returning to our numerical results, we begin with SU{12) and show in Fig. 031 how the 
average plaquette difference ((mqi ^"^^02)) varies across the transition when we vary the 'spatial' 
volume, L2. (We expect the plaquette difference to be less sensitive to domain formation than 
the Polyakov loop.) It is clear that the transition does become much sharper when we pass 
from L2 = 10 to L2 = 40 although the nature of the change between L2 = 40 and L2 = 80 
is less clear. The evidence is for a would-be first-order transition inhibited by the domain 
formation described in the previous paragraph. 

Turning now to the A^-dependence of the transition, we show in Fig. 1261 how varies 
with 7 on a 2 X 4 X 10 lattice for SU{6), SU{12) and SU{24) gauge theories. We see a 
rapid sharpening of the transition with increasing which leaves little doubt that there is a 
first-order transition at A^ = 00 at L2 = 10, and presumably at other values of L2 as well. 

4.4.3 Third transition 

To search for a third transition, characterised by a non-zero expectation value for |//x=2|, we 
take our SU{12) gauge theory on an 2 x 4 x 10 lattice and increase 7 beyond the values 
associated with the transitions discussed above. In Fig. EZI we plot the resulting values of 
(1^/^=2 1)- We see a transition of the kind that we are looking for, but one which is very smooth. 
Increasing A^ to A^ = 24 we see what appears to be a significant sharpening of the transition, 
suggesting that it might become an actual phase transition at A^ = 00. 

In Fig. I2H1 we show a histogram of the values of \l^=2\ obtained in SU (24) on a 2 x 4 x 10 
lattice at 7 = 156.25. This shows a clear peak at low values that one naturally interprets 
as belonging to the confined phase, and a further peak (or peaks) at larger values that one 
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naturally associates with the deconfined phase. This suggests that if this is a phase transition 
at = cxD then it is first order. 

5 Conclusions 

We have shown that there is a very close match in the behaviour of several observables across 
the bulk transition that separates strong and weak coupling in 2+1 dimensions, and the 
Gross- Witten transition in 1+1 dimensions. In particular the third derivative of the partition 
function, C3 oc N'^d^ logZ/dp^, appears to develop a discontinuity as — >■ 00, just as it does 
across the Gross- Witten transition, providing strong evidence for a third order transition in 
the large- limit of D = 2 + 1 SU (N) gauge theories. 

When we expressed log Z/ dp^ as a cubic correlator of plaquettes, we saw that the 
discontinuity arose from fluctuations of plaquettes at the same position. It is thus a genuine 
N = 00 phase transition that arises from the A^^ — > 00 degrees of freedom on each plaquette 
rather than from the collective behaviour of a large number of separated plaquettes. This 
motivated us to study the eigenvalue spectrum of the plaquette. We found that at the critical 
inverse 'tHooft coupling, 7 = 7^, the spectrum develops a gap at the boundary of its range 

_ g^^^ ^Yiis gap grows as 7 increases. While the gap formation does not, in itself, lead 
to a nonanalyticity in Z, it possesses a feature that does. At 7 = 7c and for A^ — > cx) the 
spectrum p{a) approaches its end-points with a vanishing derivative. This means that the 
extreme eigenvalues possess fluctuations that diverge compared to the 0(1/A^) fluctuations of 
the eigenvalues in the bulk of the spectrum, and this is directly related to the singularity in 
the partition function. All these features are exactly the same as in the D = 1 + 1 Gross- 
Witten transition. In addition we find that there is a very close match in the behaviour of 
the plaquette eigenvalue density and in the ratio of plaquette eigenvalue fluctuations when we 
compare the transitions in D = 2 + 1 and in D = 1 + 1. Thus it would appear that the bulk 
transition in 2+1 dimensions is very much like the Gross-Witten transition. 

However, there is clearly more than this going on. The Gross-Witten transition has no 
peak in the specific heat, but we see in 2+1 dimensions a clear peak that coincides with (or 
is very close to) the third order transition. The contribution from neighbouring or nearly 
neighbouring plaquettes appears to grow with A^, indicating a possible second-order phase 
transition in the large- A^ limit. Whether this is due to a correlation length that diverges as 
A^ 00 (we see a slight decrease in the lightest mass that couples to the plaquette when we 
go from SU{6) to SU{12)) or to the plaquette fluctuations decreasing more slowly than 1/A^^ 
at the critical point, is not clear at present. In any case, the fact that the correlations between 
nearby plaquettes behave as if due to a flux emerging from an elementary cube, suggests 
that the transition may be due to centre monopole(-instanton) and vortex condensation. It is 
therefore plausible that this (possible) second-order phase transition is connected to the line 
of specific heat peaks in the fundamental-adjoint plane found in SU{2) [22], which may also 
become a line of second-order phase transitions in the large-A^ limit, and which, just as in 
D = 3 + 1 can be understood in terms of condensation of Zjsf monopoles and vortices. In 
D = 3 + 1 this phase structure is believed to lead to the observed first order bulk transition. 
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Our N = oo second-order transition would appear to be a manifestation of the same dynamics, 
but in one lower dimension. Numerical calculations that are both more accurate and extend 
to larger N are clearly needed here. 

We have also investigated the sequence of finite volume transitions on L0L1L2 lattices. We 
argued that when the tori are strongly ordered, Lq <^ Li L2, these can be understood in 
terms of deconfinement, followed by high-T dimensional reduction as /3 is increased, followed 
by deconfinement in the reduced system, and so on. So the first transition, which is first-order 
for large N, occurs at /? = Pco where a{Pco)Lo = l/^co- Then as we increase /?, and hence T, the 
system will eventually be dimensionally reduced, L0L1L2 L1L2 for T Tc^^. This Li <^ L2 
system will undergo a deconfining transition when (3 = Pci Pco where a(/5c jLi = l/Tc^. 
Due to the fragmentation of the high-T phase into domains (a feature of 1 spatial dimension) 
this transition is a finite cross-over at finite N. As — > cxo the domain 'wall' tension should 
diverge (probably as N^), the domain structure will be suppressed, and the cross-over appears 
to become a genuine first order transition. At higher j3 and hence higher j'D=i+i^ ^^^j^ again 
expect dimensional reduction to occur L1L2 — > L2 at some IS some 

evidence that the 'infinitesimal' Lq x Li system undergoes a transition at a{j3c2)L2 = l/Tcj 
with > Pci- 

The first of these finite volume transitions becomes a phase transition as Li,L2 00 
at fixed A^. It also becomes a phase transition at fixed Li,L2 a.s N ^ 00. The latter will 
occur even as Li,L2 —>■ Lq. The second transition, which is a deconfining transition in the 
D = 1 + 1 effective theory, is a sharp crossover at finite A^ and appears to be first-order. As we 
remarked above, an actual phase transition is not possible at finite A^ in 1 spatial dimension 
because of domain formation. However, as we pointed out, this does not preclude a first-order 
phase transition at A^ = 00, as suggested by our numerical computations. Such a A^ = 00 
phase transition will continue to occur for L2 Li, and it may be that it also survives the 
Li,L2 —>■ Lq limit, although we have not investigated this possibility. Finally we saw some 
numerical evidence for a third N = 00 transition in the effective D = + 1 theory, although 
here the calculations are no more than suggestive. These arguments can trivially be lifted 
to SU{N) gauge theories in 3+1 dimensions, where they clearly have some relation to the 
N = 00 finite volume transitions on lattices transitions discussed in [21 El El • 

In view of conjectures [21 IHIIZI that Wilson loops in D = 3 + 1 may undergo N = 00 non- 
analyticities, when their area, in physical units, reaches a critical value, we have analysed the 
behaviour of Wilson loops in Z) = 2 + 1 SU{N) gauge theories. Our results show a remarkable 
match between the behaviour of Wilson loops in D = 2 + 1 and in D = 1 + 1, where a gap 
in the Wilson loop eigenvalue spectrum is known to open at a critical area at A^ = 00, in 
both the lattice and continuum theories [HI GDI • We find that the eigenvalue spectra of Wilson 
loops (and indeed Polyakov loops) in D = 2 -|- 1 match those of Wilson loops in D = 1 + 1 
when the traces are equal. Moreover the spectra of Wilson loops of any size (in lattice units 
and when larger than about 2x2) also match if the couplings are tuned to values where their 
traces are equal. This is true for any fixed A^. As a corollary, it immediately follows that in 
D = 2-|-latA^ = ooa gap will form in the eigenvalue spectrum of a Wilson loop at a critical 
coupling that depends on the size of the loop. However because of a logarithmically divergent 
self-energy piece, this non-analyticity in the spectrum will not occur at a finite non-zero value 
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of the area in the continuum hmit. This is in contrast to the case in D = 1 + 1. We have 
prehminary evidence |^ for a similar matching between Wilson loops in D = 3 + 1 and those 
with the same trace in lower dimensions. Here the self-energy divergence is even more severe 
and the gap forms deep in the ultraviolet. As in D = 2 + 1 one can imagine regularising this 
self-energy by using finite mass sources in constructing the 'Wilson/Polyakov loops', so as to 
obtain a gap formation at a fixed physical area. 

The appearance of the gap at iV = oo follows quite generally if we make plausible con- 
nections with Random Matrix Theory. The spectrum of an / x / = aL x aL Wilson loop 
should be flat at large a, since linear confinement demands its trace to be oc exp{— cr/^} ~ 0, 
while at sufficiently small a we expect to find the Wigner semi-circle of the N = oo Guassian 
Unitary Ensemble. Somewhere in between a gap must form. Because the derivative of the 
spectrum diverges at its end-point, in contrast to the plaquette at the bulk transition, there 
are no anomalous fiuctations of the extreme eigenvalues and no non-analytic behaviour in the 
correlators that are related to derivatives of the Wilson loop with respect to the coupling. 
And indeed we find the traces of Wilson loops to be analytic in the coupling just as they are 
in D = 1 + 1. Thus the physical implications of this nonanalyticity in the eigenvalue spectrum 
remain unclear. 

The remarkable similarity between the eigenvalue spectra of Wilson loops in different 
dimensions does not appear to have a simple explanation within Random Matrix Theory and 
merits a more careful and quantitative investigation than the one provided in this paper. 

The D = 2 + 1 large- phase structure that we have investigated in this paper can be 
understood, as we have argued above, in terms that appear to allow a unified understanding 
of these phase transitions mD = l + l, D = 2 + l and D = 3 + 1 SU{N) gauge theories. 

Note added. 

This revised version arose from our discovery, immediately after sending the original version 
to the archive, of the papers HH] which then motivated our revised and extended study of 
Wilson loops in this paper. As this revision was in progress an interesting paper [20] on gap 
formation in smeared Wilson loops in D = 3 + 1 has appeared. 
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N 


7c 


2 


0.700(3) 


6 


0.719(3) 


12 


0.722(2) 


24 


0.722(1) 


48 


0.722(2) 



Table 1: Inverse coupling at which gap forms for 3 x 3 loops in SU{N). 




Figure 1: The spectrum of eigenvalues, e*", of a 3 x 3 Wilson loop for SU{48) in D = 1 + 1 
at the critical coupling 7 = 1/A = 1.255, together with the continuum spectrum ( ). 
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Figure 2: The average plaquette as a function of 7 = = 2^ for SU(6) (+), SU(12) (x), 
SU(24) (*) and SU(48) (□). 
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Figure 3: The specific heat, C2, as a function of 7 = — ^ = for SU(6) (sohd fine) and 
SU(12) (dashed hne). 



30 




0.6 0.65 



Figure 4: The 'local' specific heat, P2, as a function of 7 = 3^ for SU(6) (+), SU(12) (x), 
SU(24) (*) and SU(48) (□). 
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Figure 5: The specific heat, C2 (equal to P2 here) as a function of 7 = in 1+1 dimensions 
for SU(6) (+), SU(12) (x), SU(24) (*) and the analytic result for SU(oo) (solid line). 
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Figure 6: The cubic local plaquette correlator, P3, as a function of 7 = for SU(6) (+), 
SU(12) (x), SU(24) (*) and SU(48) (□). 




Figure 7: The cubic plaquette correlator, C3 (equal to P3 here) as a function of 7 = 2^ in 
1+1 dimensions for SU(6) (long dashes), SU(12) (short dashes), for SU(6) (+), SU(12) (x), 
SU(24) (*) and the analytic result for SU(oo) (solid hne). 
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Figure 8: Fluctuations of extreme plaquette eigenvalues, i?p, as a function of 7 
SU(6) (+), SU(12) (x), SU(24) (*) and SU(48) (□). 
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Figure 9: As in Fig. |Hl but in 1+1 dimensions. 
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Fi gure 10: Density of plaquette eigenvalues, e'", for SU(12) in 1+1 dimensions at 7 — — 
0.462 (long dashes) and 7 = 0.542 (dots) and in 2+1 dimensions at 7 = 0.417 (solid line) and 
7 0.451 (short dashes). 




Figure 11: The plaquette correlator, Co, as a function of 7 = for SU(6) (+), SU(12) (x), 
SU(24) (*) and SU(48) (□). 
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Figure 13: Trace of 2 x 2 Wilson loop in SU(6) (+) and SU(12) (x) and of the 3 x 3 loop in 
SU(6) (★) and SU(12) (□). 
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Figure 14: Local 'specific heat' of 2 x 2 and 3x3 Wilson loops in SU(6) and SU(12). Labels 
as in Fig. [T21 
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Figure 15: Quadratic correlator of 2 x 2 Wilson loops, W2^^, as a function of 7 = for 
SU(6) (+), SU(12) (x), SU(24) (*) and SU(48) (□). 
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Figure 16: Cubic correlator of 2 x 2 Wilson loops, W^^ , as a function of 7 = ior SU(6) 
(solid line) and SU(12) (long dashes). 




Figure 17: i?^^^ as a function of 7 = ^ for SU(6) (solid line), SU(12) (long dashes) and 
SU(24) (short dashes). 
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Figure 18: 3x3 Wilson loop eigenvalue density, e*", for SU(12) in 1+1 dimensions at 7 = 
— 1.255 (solid line) and in 2+1 dimensions at 7 = 0.722 (long dashes), and the continuum 



2N 



large-N distribution in 1+1 dimensions aX A = Acru (short dashes). 




Figure 19: 3x3 Wilson loop eigenvalue density, e*", for SU(12) in 1+1 dimensions at 7 = 
2^ — 2.215 (sohd hne) and in 2+1 dimensions at 7 = 1.111 (long dashes), and the continuum 
large-N distribution in 1+1 dimensions a,t A — 0.539 Acru (short dashes). 
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Figure 20: Eigenvalue density in SU(6) in 2+1 dimensions for the 2x2 loop at 7 = 0.483 
(solid line), the 3x3 loop at 7 = 0.719 (long dashes), the 4x4 loop at 7 = 0.965 (short 
dashes) and the continuum large-N distribution in 1+1 dimensions at A — Acru (dots). 




Figure 21: Couplings for which the gap forms for L x L Wilson loops in SU(2) (+) and SU(6) 
(*), and fit to SU(2) data (dashed line). 
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Figure 22: Polyakov loop eigenvalue density in SU(12) in 2+1 dimensions in the confined 
phase at 7 = 0.764 (solid line) and in the dcconfined phase at 7 = 0.833 (long dashes), and 
the 3x3 Wilson loop in 1+1 dimensions at 7 = 1.684 (dashes). 
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Figure 23: Values of the shorter 'spatial' Polyakov loop (|/^=i|), •, and the plaquette difi^erence 
500 X {{U01 — U02)), o, on a 2 x4x40 lattice in SU(12) versus the bare inverse 't Hooft coupling, 
7 = P/2N'^. 
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Figure 24: Histogram of values of the spatial Polyakov loop, \lij,=i\, in SU(12) on a 2 x 4 x 80 
lattice at 7 = 3.368. 
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Figure 25: The average plaquette difference {6u) = 10'^ x ((mqi ~ ''^02)) in SU(12) on 2 x 4 x L2 
lattices with L2 = 10 (+), L2 = 40 (o), and L2 = 80 (•), versus 7 = p/2N^. 
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Figure 26: The average n - 1 Polyakov loop for SU{<o) (+), SU {12) (o), and 5?7(24) 
versus the inverse bare 't Hooft couphng 7 = I3/2N'^, all on 2 x 4 x 10 lattices. 
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Figure 27: The average fi = 2 Polyakov loop in SU{12), o, and in 5'C/(24), •, on a 2 x 4 x 
lattice versus the inverse bare 't Hooft coupling 7 = I3/2N'^. 
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Figure 28: Histogram of values of \li^=2\ in SU(24) on a 2 x 4 x 10 lattice at 7 = 156.25. 
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